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We calculate in perturbative QCD the transverse momentum imbalance of dijet and dihadron 
production in high energy p+A (d+A) collisions. We evaluate the effect of both initial- and final- 
state multiple scattering, which determines the strength of this transverse momentum imbalance. 
Combining this new result with the suppression of the cross section in d+Au collisions, which arises 
from cold nuclear matter energy loss and coherent power corrections, we are able to describe the 
dihadron correlations measured by both PHENIX and STAR collaborations at RHIC, including 
c I ^ , mid-mid, mid-forward, and forward-forward rapidity hadron pairs. 
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PACS numbers: 12.38.Bx, 12. 39. St, 24.85. +p, 25.75.Bh 



, In recent years, there has been renewed interest in nuclear effects that affect hadronic observables in proton- 
O i' nucleus (p+A) and deuteron- nucleus (d-l-A) collisions [ij. Much of this interest was sparked by forward rapidity 
measurements in d-|-Au collisions at the Relativistic Heavy Ion Collider (RHIC) that indicate considerable suppression 
of both single hadron and dihadron production cross sections relative to the binary collision scaled proton-proton 
(— I ' (p+p) baseline Back-to-back dihadron correlation measurements have attracted the most attention and are 

'— ' : under active investigation in the framework of different theoretical formalisms 04111 ■ On one hand, quantifying the 
differences between p-|- A (d-|-A) and p+p collisions (referred to as "cold nuclear matter effects" ) can provide a solid 
baseline for the unambiguous identification of any additional final-state hot and dense medium effects in heavy ion 
collisions [HI . On the other hand, the study of p+A (d+A) collisions at forward rapidity is interesting in its own right 
^sj since it probes the inelastic and coherent multiple parton scattering in the nucleus, as well as any possible modification 
■ of the small-x parton distribution functions }7Hlll [T^-EB] ■ 

I In perturbative QCD, at leading order, high transverse momentum {p±) particle production arises from hard 2 — >■ 2 
. scattering processes. In this picture, two incoming collinear partons scatter on each other to produce two high 
transverse momentum partons, which then fragment into jets of hadrons. These jets, as well as the hadron pair 
formed by their leading particles, are approximately back-to-back in the transverse plane. In high energy proton- 
nucleus reactions the incoming partons that participate in the collision can undergo multiple interactions. Even in 
the presence of multiple parton scattering, particle production at high transverse momentum is dominated by a single 
hard interaction (x p]_ ^ where is the typical scale of the soft transverse momentum exchanges. We refer 
to the interactions that precede the large scattering as initial-state and the interactions that follow the large 



. scattering as final-state. 

Both initial- and final-state multiple interactions can affect dijet (dihadron) production in p+A (d+A) reactions. 
More specifically, they lead to an increase in the transverse momentum imbalance jlTj] that can be perturbatively 
computed within a high- twist formalism [iGl - fTsj . In this paper we give a theoretical derivation of the increase in 
the transverse momentum imbalance and demonstrate that it will lead to a broader away-side peak for the dihadron 
azimuthal correlation distribution in d+A collisions, consistent with the experimental observations at RHIC. Fur- 
thermore, multiple parton scattering also manifests itself trough medium-induced radiative corrections [l^, [2^ and 
power suppressed contribution to the cross sections. In reactions with nuclei (p+A and A+A) initial-state radiative 
energy loss effects are always present but their effect is most pronounced at forward rapidity [l^ HIl ■ In the coherent 
scattering regime, the high-twist contributions are enhanced by the nuclear-size cx A^^^ and may become important 
at small and moderate transverse momenta. Such nuclear-size enhanced power corrections have been resummed for 
both inclusive single hadron and dihadron production processes in p+A collisions ■ In the phenomenological part of 
this paper we combine the broadening in the away-side width with the nuclear cross section suppression to study the 
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modification of the production rate and observed shape of dihadron azimuthal correlations in going from p+p to p+A 
(d+A) collisions. With parton scattering parameters in cold nuclear matter constrained by deep inelastic scattering 
experiments [13] , calculations are consistent with the experimental measurements at RHIC 0| . 

The rest of our paper is organized as follows: in Sec. II we take into account both initial- and final-state multiple 
parton scattering to calculate the increase in the transverse momentum imbalance of dijet and dihadron production 
in p+A (d+A) collisions. In Sec. Ill we first evaluate the width of the away-side peak of the dihadron correlation 
distribution at RHIC and show that our formalism can describe the PHENIX experimental data well if we take 
into account the nuclear-induced increase in the transverse momentum imbalance. We then overview the nuclear 
effects that lead to the suppression of single and double inclusive hadron production in p+A (d+A) reactions. We 
demonstrate that the calculated mid-forward and forward-forward rapidity hadron pair production attenuation is 
consistent with recent PHENIX measurements. At the end of this section, by combining the nuclear suppression 
and the increase in the transverse momentum imbalance, we are able to describe the dihadron azimuthal correlation 
distribution observed by the STAR experiment at RHIC. Summary and conclusions are presented in Sec. IV. 



II. DIJET AND DIHADRON TRANSVERSE MOMENTUM IMBALANCE IN p + A COLLISIONS 

In this section we derive the increase in the transverse momentum imbalance of back-to-back jet production in p+A 
collisions. This phenomenon is often referred to as nuclear-induced broadening in the literature. We then generalize 
the formalism to study back-to-back hadron production. 



A. Dijet transverse momentum imbalance 



We start by describing the pQCD formalism for evaluating dijet production in p+p collisions: 

p{P') +p{P) ^ Jl(Pli) + J2(P2±) + A. (1) 

Let us define P± = \Pi± — P2±\/2 to be the magnitude of the average transverse momentum of the jet pair. In leading 
order perturbative QCD, the jets are produced back-to-back, Pij^ = —P2±, thus \Pi±\ — \P2±\ = P±- The differential 
cross section can be written as 1231: 



da 



dyidy2dP 



_ ^al ^ fa/p{x')h/p{x) 

a,b 



XX 



(2) 



where X)a b runs over all parton fiavors, s ~ {P' + P)^, fa,b/p{x) are the parton distribution functions with the 



(e ^1 + e y^), and yi, y2 the rapidites of the two jets. H^t,_^cdi^' 



momentum fractions x' — (e^^ + e^^) , 

are the partonic cross sections as a function of the usual partonic Mandelstam variables s,t,u. These cross sections 
are well-known [23|. we will reproduce them here for later convenience: 
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where TVc = 3 is the number of colors. 




FIG. 1. Left: Feynman diagram for the partonic qa[ — >■ qq scattering channel. Right: initial-state double scattering Feynman 
diagram that contributes to the dijet transverse momentum imbalance increase. 



Let us now study dijet production in proton-nucleus collisions. In such collisions the energetic incoming parton from 
the proton can undergo multiple scattering with the soft partons inside the nuclear matter before the hard collisions. 
We refer to these interactions as initial-state multiple scattering. After the hard collisions, the two leading partons, 
produced almost back-to-back, will eventually hadronize into jets as observed in the experiments. However, these 
partons will also likely undergo multiple interactions in the large nucleus. We refer to these interactions as final-state 
multiple scattering. Such separation is only possible when the scale of hard scattering is considerably larger than the 



scale of typical transverse momentum exchanges between the projectile and the medium, ^ ^ [2l[. The initial- 
and final-state multiple scattering lead to dijet acoplanarity, or imbalance. To quantify this effect, let us define the 
dijet transverse momentum imbalance g_L as: 

gl = Pi± + P2j., (11) 

and the averaged transverse momentum imbalance square as: 

(gi) = / rf^glgi , ■ I , (12) 

J dyxdyidPl^d^q^l dyidy2dPl 

The nuclear enhancement (or broadening) of the transverse momentum imbalance in p-f A collisions relative to p+p 
collisions can be quantified by the difference: 

A(gi) - {ql)pA - {ql)pp. (13) 

As demonstrated in Ref. [161] . the transverse momentum imbalance increase A(g^) can be calculated perturbatively 
and the leading contribution comes from double scattering. To illustrate the method, we will study the simple partonic 
channel qq' — > qq' . At lowest order, the partonic cross section H^^,^^^, is calculated from the Feynman diagram in 
Fig- Q] (left). At this order, the transverse momenta of the two outgoing partons (or the dijet) are equal and opposite 
and we have: 

'^'^ oc S\q^). (14) 



dyidy2dPld'^qi_ 



The leading contribution to the transverse momentum imbalance increase in p-f A collisions comes from either initial- 
state double scattering, as in Fig. [1] (right); or final-state double scattering, as in Fig. [51 Let us first focus on the 
calculation of initial-state double scattering. In this case, the small fcj_-kick in the nucleus will generate a small 
transverse momentum imbalance for the dijet: 

=Pi±+P2± =fc±. (15) 
Following [l^, we evaluate the contribution from the double scattering diagram to {q\)pA'- 

/■ ,2- 2 dapA T^alfq/pix') /■ ,2- 2 /■ , , ,2, 7p / , ^ 

j d q^q^ dj; dy dP^d'^q ^ — J '^^'^^ J '^^I'^^^a k±_T Aq{x,Xi,X2,kA^) 

xH{x, Xi,X2,k±, x'p')6'^{q± — k±) , (16) 
where the matrix element TAq is represented by the bottom blob in Fig. [T] (right) and has the following expression: 

Ta (t Xo k,\— f '^^ f !^^y^ pixip'^Vi pi(x-xi}p+y- -i(x-X2)p+y^ -ik±-y± 

iAg[x,xi,X2,ii±)- J 2^ J ^2^^2e e e e 

X 1 {pa\A+ {y^ , 0^)7^(0)7+ V(j/i" )A+ {y- , y±)\pA) ■ (17) 
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As in [i^ the calculation first takes advantage of 5'^{q± — k±) to set = fc^. These two factors of k± can be 
converted into transverse derivatives on the two fields in matrix element Taq in Eq. ([T7| by performing a partial 
integration for cfy±: 

klA+{y^,0^)A+{y-,y^)^F+{y^,Q^)F+"{y-,y^). (18) 
We then expand in k± in the partonic part H around fc^ ~ 0, and keep the first non-vanishing term: 

xi, X2, xV) = ^(2;, xi, 3:2, fcj^ = 0, x'p'). (19) 

In other words, since the fields have been converted to the gauge-covariant gluon field strength in the matrix 
element T^q, one can set fcj^ = in the hard part. One will later find that all such hard-parts H{x,xi,X2,x'p') for 
the double scattering diagrams reduce to the lowest order Born diagrams up to color factors. Thus, the calculation of 
transverse momentum imbalance enhancement A((7^) is simpler than the usual higher- twist calculations in which one 
calculates the multiple scattering contribution to the differential cross section itself directly (instead of the gr^-weighted 
one). 

We obtain: 

2^ 2 (^^pA nal fq/pix') , 
^'''^-'^ dy dy dP^cPq± ^ a/x — J dxidx2TFq[x, xi, X2)H(x, xi, X2, x p ), (20) 

where Tpg is a twist-4 four-parton correlation function: 



t(;)(x,xi,X2) ^^e^^^P^y^e^^^--^^P^y-e-^^^-^^^P^y^^{pA\F^+{y^-)^^^^^ 

Here and thereafter we will use the superscript "I" ("F") to indicate the contribution associated with initial- (final-) 
state multiple scattering. The hard partonic function i/(x, xi, X2, x'p') is given by: 



H{x,xi,X2,x'p') = Stt a 
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27r xi — X — X2 — X -|- 



where Cp ~ {N^ ~ i)/2Nc. Substituting Eq. ((22)) into Eq. pop and performing the integration over xi and X2, we 
obtain: 

dyi dy2 dPl d^q± ~ 1 J 



where r^^^]^(x) is the twist-4 quark-gluon correlation function defined as (TgI. [iTj: 

^,m(^) = / ^e""^'"/ ^^^^(y"-yr)^(-y2")^(pA|F„+(y2-)Vi,(0)7+V',(y-)^^+"(2/r)bA), (24) 

and Hqq'^qqi is the hard-part function given by: 

Hqq'^qq' ^ ^pHqq'^qq'- (25) 



±2^ 



FIG. 2. Final-state double scattering Feynman diagrams that contribute to the dijet transverse momentum imbalance increase, 
for the partonic channel qq — > qq' . 



Next, we calculate the contribution to the transverse momentum imbalance from the final-state double scattering 
diagrams in Fig. [21 This evaluation is slightly more complicated since there are four diagrams as opposed to a single 
diagram for the initial-state case. However, the calculation is very similar and straightforward. Different Feynman 



5 



diagrams only differ in terms of color factors and their sum provides a measure of the average square transverse field 
strength probed by the outgoing partons .17]. We obtain: 
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( Sn^as \ Tra^ fq/p{x') >f) 



dy^dy^dPlcPq^ ~ I J s 



(26) 



where T^J^ is given by the same expressions in Eq. (p4|) . except for the 6'- functions that are replaced as follows [T6l.[T7|: 

e{y- - y^)9{-y^) ^ e{y^ - y-)9{y^) 
and the final-state hard-part function is given by 



(27) 
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2m 



(28) 



Likewise, we can go ahead to calculate the contribution to the transverse momentum imbalance from both initial- 
and final-state double scattering diagrams for all partonic channels. The calculation is straightforward though tedious, 
the result can be summarized as: 



72 ^ 2 

d qj_q_i^ 



da^ 



pA 



dyidy2dPld^qi_ - I J s 



2 ^^^^^ rftf rp 



TlpA^)Hi,^,^is, t, u) + T^j;{x)H:,^,^is, i, u)\ ,(29) 



where t'^^^'^\x) are given in Eqs. (|24p and (P7|. In the calculation for the partonic channels qg — > qg and gg — > gg, 
see for example the Feynman diagrams shown in Fig. [31 two other similar twist-4 gluon-gluon correlation functions 
Tj^^^'^\x) appear. The operator definition of r^'^]^(x) is given by (l6l. [T7|: 



while TgJ^{x) is given by the same expression with the 0- function replacement specified in Eq. ([27l 
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FIG. 3. Sample diagrams that contribute to the dijet transverse momentum imbalance for the partonic channels qg qg (left 
two) and gg — !> gg (right two), where (a) and (c) are for initial-state double scattering, (b) and (d) are for final-state double 
scattering. 

Finally, from Eqs. ([2]) and ([29]), we obtain the transverse momentum imbalance increase A{q'j_) in pA collisions: 

Ea,^ \TiPAx)Hi,^^,is, i u) + TiJl[x)H^,^^,{s, i, u) 



/a/p(^')/Vp(a:) ttV I- 1 f,\ 
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(31) 



The partonic hard-part functions Hj^i^^^^ are associated with the initial-state multiple scattering, and are given by: 
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CAH^b^^d a=gluon 

with Ca = On the other hand, H^f^^^^ are associated with the final-state multiple scattering, and are given by 
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Dihadron transverse momentum imbalance 



One can easily generalize the transverse momentum imbalance for dijet production to dihadron production. To start, 
we write down the leading order differential cross section for dihadron production in p+p collisions, p{P') + p{P) — > 
hi{pi±) + h2{p2±) + X, 



da 



dyidy2dpij_dp2j_ 



2tU^^ [ dzi Ja/p{x')fb/p{x) jj f 



ibcd 



where Df^^/^{zi) and Z3/j2/d(^2) are fragmentation hmctions, Z2 = ZiP2±/pi± p^ . and: 



Pi± 



Pi± 

Zly/s 
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(41) 



(42) 



At this order, the final hadron pair comes from the fragmentation of the back-to-back parton pair in the partonic 
collision ab — cd. One has pi± — ziP± and p2± — —Z2P± with P± the transverse momentum of the first parent 
parton. Then, the dihadron transverse momentum imbalance q± can be written as: 



q± = Pl± +P2± = {zi - Z2)P±. 

Thus, the averaged transverse momentum imbalance {q'j_) in p+p collisions is: 

da I da 



J dyidy2dpi 



±dp2±d'^q±/ dyidy2dpi±dp2± 



{{z,-Z2fPl). 



(43) 



(44) 



On the other hand, in p-l-A collisions there will be transverse momentum imbalance increase generated by the small 
A:_L-kick from the initial- and final-state multiple scattering in the nuclear medium, as demonstrated in the previous 
section. Since the two effects are independent: 



{<ll)pA = {{{Zl - Z2)Pl_ + kl_f) = ((Zl - Z2fPl) + (fci). 



(45) 



where in the second step we have taken ((zi — Z2)P± ■ k±) — due to the random nature of the fc_L-kick. Thus, the 
transverse momentum imbalance increase in p-l-A collisions when compared to p-f p collisions is given by: 



A{ql) = {ql)pA - {ql)pp = (kl), 



(46) 



and is equal to the strength of the soft fcj^-kick from the multiple scattering. From the result of the last section, we 
immediately obtain for the dihadron production: 
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(47) 



Eqs. pil) and (|47l) are the main new theoretical results of our paper. We will use them in the phenomenological 
studies of dihadron azimuthal correlations in the next section. 
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III. DIHADRON CORRELATION IN d + Au COLLISIONS 



In this section, we study the phenomenological apphcations of our resuhs. First, we use the transverse momentum 
imbalance result to derive the width of the away-side peak in dihadron correlations and compare our findings to 
existing experimental data. Next, we review the nuclear effects which lead to the suppression of dihadron production 
in d-|- Au collisions and calculate the nuclear modification factor w^^J^ , which turns out to describe the data reasonably 
well. Finally combining the calculation on both the away-side width and the nuclear suppression factor, we are able 
to describe the broadening of the dihadron azimuthal correlations in d-|-Au collisions relative to p-l-p collisions. 



A. Dihadron imbalance and the width of the away-side peak 

In order to estimate the nuclear- induced broadening for dihadron production given by Eq. (|T7| , one needs to know 
the four-parton correlation functions Tj^j^^x) with b — q,g. Following Refs. 0, H^l, these can be modeled as: 



Such decomposition into a leading twist parton distribution function and nuclear enhanced transverse momentum 
transfers squared and/or power corrections is motivated by the 0-function structure, and the similarity of the operator 
definition between T^^^JJ^{x) and fq^g/Ai^)- In Eq. pS)) = 0.09 — 0.12 GeV^ represents a characteristic scale of 
high twist corrections per nucleon and has been extracted from deep inelastic scattering data [22^ It has also been 
employed to describe single inclusive hadron production in d-|-Au collisions at forward rapidities 0, [l^ [2l| . We can 
now write the increase of transverse momentum imbalance as: 

The transverse momentum imbalance {q'^)dA — {q±)pp + characterizes the shape of the azimuthal correlation 

of back-to-back dihadron production in d-|-Au collisions, as it is closely related to the width of the away-side peak 
ap- First, we have the following relation [2^ [25|: 



(bouti)' = {\j±y\r+xl {{\j±y\r+2{\k^y\)l,^,J , (50) 

where (|pout|) is the average transverse momentum out of the plane defined by the momentum of the trigger particle 
RL.tiig and the beam axis, j±y is the component of the particle momentum perpendicular to the jet momentum in 
the fragmentation process, and (|fc_Lj,|) parton = (^1) parton /tJ". On the other hand, partonic (fc^)parton is related to the 
transverse momentum imbalance as {q'j_) — 2(fc^)parton. Realizing that: 

(bouti) =P_L,assoc sin (|A(?!)|) =p_L,assoc sin ^y'^CTF^ , (51) 
where pj^.assoc is the transverse momentum of the associated particle, and 

7^_L,assoc ' P_L.triff P_L.assoc /i * i\\ J'_L,assoc \ \ " \ ( r c\\ 

XE = — ~ ■ COS (|A^|) = ■ cos I \j -ap ) , (52) 

Pj_,triK PJ-.trig -P-L,trig 




we obtain: 




{\3±v\r 



— 1 — ~r 5 

^_L,trig ^_L,trig 



(53) 



It has been shown that ~ 400 MeV is sensitive only to vacuum fragmentation and is independent of the 

center-of-mass energy ^/s and the trigger particle momentum 0, HI]. Thus, for the selected trigger and associated 
particles (with specific momenta), the width of the away-side peak ap depends closely on the transverse momentum 
imbalance {q\}. 

The PHENIX collaboration at RHIC has measured ap, as shown in Fig. |31 We first compute A(g^) by taking 
= 0.12 GeV^ [7] and by using the CTEQ6L1 parton distribution functions [26] and the deFlorian-Sassot-Stratmann 
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J \ I I I \ \ \ \ I \ \ \ L 



1 2 3 4 5 

PT,a.soc(GeV) 

FIG. 4. The width ctf of the away-side correlation is plotted as a function of the associated hadron transverse momentum 
p±, assoc. Both hadrons are at mid-rapidity \yi.2\ < 0.35. The upper panel is for 3 < px,trig < 5 GeV and the lower panel is 
for 5 < px,trig < 10 GeV. Solid curves represent the width calculated from {qj_)dA ~ {<1j_)pp + ^(^l) a-nd include multiple 
scattering effect. Dashed are calculated from {q±)pp. Data is from PHENIX j^. 



(DSS) hadron fragmentation functions [l^- We choose {q'±}pp ~ 3.3 GeV^ extracted from PHENIX p-l-p data [13] to 
evaluate {q\)dA in d-|-Au collisions. 

In Fig.mthe width ap of the away-side peak is plotted as a function of the associated particle transverse momentum 
P_L,assoc for two ranges of the trigger particle momentum in the minimum bias d-|-Au collisions. The solid curves are 
calculated from {q\)dA, while the dashed ones are from {q\)pp (i.e., without multiple scattering). Even though the 
increase in the transverse momentum imbalance ^{q\) is on the order of 2 — 3 GeV^, the broadening effect in the width 
(Ti? is moderate for the mid- mid correlation, where both trigger and associated particles are in the mid-rapidity region 
\yi.2\ < 0.35. Nevertheless, the calculated ap with multiple scattering effects appears to agree with the PHENIX data 
slightly better. 



B. Nuclear modification factor in d+Au collisions 



Coherent multiple parton scattering can also affect the rate of dihadron production. The difference between p-l-A 
(d-|-A) and p+p reactions is usually quantified by the nuclear modification factor: 

^(2) _ dadA/ dyidy2dpij_dp2± ^^^^ 
'^^ {Ncoii)dapp/dyidy2dpi_Ldp2± 

There are two major effects that control the magnitude of -R^^: dynamical shadowing [T3| and cold nuclear matter 
energy loss [l^, |2l|, [2^. Dynamical shadowing effects have been shown in Refs. Il3| to contribute to the cross 
section at the power corrections level. Nuclear size enhanced (A^^^) power corrections can be resummed for a given 
partonic channel, as shown by Qiu and Vitev, and lead to the following shift in the momentum fraction x for the 
parton inside the nucleus (t-channel): 

X X (^l + Cd^^^-^l^y (55) 

where Cd — ^ (Ca/Cf) if parton c? is a quark (gluon) in the partonic channel ab — >■ cd. The reason for the values 
of Cd is that has been determined for coherent quark scattering [22| . For the hard scattering partons that couple 
to the nucleus through other channels (c in the u-channel and a in the s-channel), similar shifts in the momentum 
fraction x have been derived in These can be obtained easily from Eq. (I55p by substituting t u and t s, 
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respectively, and by keeping track of the flavor of partons c and a: Cd — Cc and Cd — >■ Ca- In our paper we have 
taken into account the shift in x consistently for all three channels. 

It is important to emphasize the similarities and differences between the dynamical shadowing calculation performed 
by Qiu and Vitev and the nuclear enhancement in the dijet (dihadron) imbalance presented in our last section. Qiu 
and Vitev studied the multiple scattering contributions to the dihadron differential cross section directly. These 
contributions are power suppressed by the hard scale, for example cx A^/^S^"^ /i for the i-channel. In Refs. [3, [13] the 
coherent multiple scatterings have been resummed and have been shown to lead to the shift in x specified in Eq. (j55p . 
When this shift is small, only the first term in the resummed series is important and it is directly proportional to 
the twist-4 correlation function T^^^J^{x). On the other hand, here we calculate the multiple scattering contribution 

to the g^^-weighted differential cross section. In the evaluation of the double scattering contribution to this weighted 
cross section we can set fcj, = in the hard parts. Because of the same partonic scattering origin, the dijet imbalance 
enhancement A(g^) and the nuclear modification factor i?^^ depend on the same characteristic scale 

As the parton from the proton undergoes multiple scattering in the nucleus before the hard collisions, it can 
lose energy due to medium-induced giuon bremsstrahlung. The spectrum of this initial-state energy loss was first 
derived in [l9j . At collider energies, for mid and forward rapidities {y > 0) even particles of small p±_ come from 
partons of very high energy in the rest frame of the nucleus (A is at large negative y). In this regime, initial- 
state cold nuclear matter energy loss can noticeably affect the experimentally measured cross sections [131 ■ There is 
renewed interest in constraining its magnitude through measurements with electromagnetic final states, such as the 
Drell-Yan production [2l|, [2^ . Initial-state energy loss has been implemented in the evaluation of hadronic and jet 
observables [H, [1^, H^l • Such studies build upon early cold nuclear matter energy loss phenomenology [sTI - is^ . 

For a generic differential cross section, medium-induced radiative corrections factorize from the hard scattering 
process and enter as an integral convolution. For initial-state energy loss, a simple change of variables presents the 
effect as a rescaling of the momentum fraction of the incoming parton from the proton |,21i] : 

/(.^) ^^{\M{x'P',xPm2^)H^{p,~pf) , (56) 

where {\M{x' P' ,xP)\^) = H^b^^j^{s,i,u) defined in Sec. II. For a single (ug = 1) emitted gluon dN9{e)/de = P(e) 
is the probability distribution for fractional energy loss e. In general, ng ^ 1 and the probability distributions 
Pqg{e) is constricted from dN3{e)/de in the independent Poisson approximation, e = '^^uJi/E. In the soft gluon 
approximation medium- induced radiative corrections always factorize from the hard short-distance scattering jl9l [35| . 
For arbitrary kinematics, this factorization is always exact for parent quarks (ssj . Just like in the vacuum, beyond the 
soft gluon a ppr oximation medium-induced splitting kernels for parent gluons factorize only for polarization-averaged 
observables [20| . In this manuscript we work in the soft medium- induced splitting approximation which allows for an 
energy loss interpretation of the cross section suppression. From Eq. ([55)) . this effect is easily implemented through: 

fU^') ^ f deP,{e)U, (^^) , /g(x') ^ dePg{e)fg (^^) . (57) 

In Fig. [51 we plot the nuclear modification factor R^J^ in d-|-Au collisions for dihadron production as a function of 

mean binary collision number (Ncoii) for both mid-forward (top panel) and forward-forward (bottom panel) correlated 

pairs. To take into account the centrality dependence, in both the transverse momentum imbalance and the dynamical 

shadowing we have replaced A^/^ — > A^^^ {N^^^^{b)) / {N^^i{bniin.hias)) Q- In the calculation of initial-state radiative 

energy loss [l^ [2l| the average length of the medium is similarly scaled. For both mid-forward and forward-forward 

cases, the associated particles are in the forward-rapidity region 3.0 < 2/2 < 3.8 with momenta 1.0 < p_L,assoc < 1-5 

GeV. For the mid-forward case the trigger particle is at mid-rapidity \yi\ < 0.35, while for the forward-forward case 

the trigger particle is at forward rapidity 3.0 < yi < 3.8. The solid curves contain both dynamical shadowing and 

energy loss effects, whereas the dashed curves contain only the dynamical shadowing effect. As we can see, the 
(2) 

calculated §i^^ ^ very good description for the mid-forward correlated pairs. The current experimental data still 
has large uncertainties, and is not able to further constrain the individual contribution of dynamic shadowing and 
the cold nuclear matter energy loss. For the forward-forward correlated pairs our formalism is still roughly consistent 
with the PHENIX measurements though the agreements get slightly worse. In particular, for the most-central 
collisions the data seems to indicate a very strong suppression. Initial-state energy loss plays an important role in 
bringing the calculation closer to the data, as it does for inclusive particle production at forward rapidities [2l|. It 
will be interesting to see whether such suppression factor can be explained within other formalisms, for example those 
based on color glass condensate [1, Q . 
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FIG. 5. Nuclear modification factor i?^'^ for back-to-back dihadron production in d+A collisions. Solid curves contain both 
dynamical shadowing and energy loss effect, whereas dashed curves contain only dynamical shadowing. Both top and bottom 
panels have an associated particle in the forward-rapidity 3.0 < 1/2 < 3.8 and 1.0 < px,assoc < 1.5 GeV region. From left to 
right there are four groups corresponding to a mean binary collision number {A'^coii)=3.2, 6.6, 10.2, and 15.1, respectively, even 
though they have been offset from these actual values for visual clarity. Top panel: trigger particle at mid-rapidity \yi \ < 0.35. 
Within each centrality selection, from left to right, the transverse momentum of the trigger particle is in the [1, 2], [2, 4], [2.5, 
3], [3, 4], and [4, 7] GeV interval. Bottom panel: trigger particle at forward rapidity 3.0 < yi < 3.8. Within each centrality 
selection, from left to right, the transverse momentum of the trigger particle is in the [1.1, 1.6], [1.6, 2.0], and [2.0, 5.0] GeV 
interval. Data is from PHENIX ^. 



C. Dihadron azimuthal correlations 



It is also useful to present directly the dihadron azimuthal correlation distribution, see for example the STAR 
measurement in Fig. [B] In this case, the dihadron correlation can usually be approximated by two Gaussians for the 
near-side and the away-side, and a constant background: 

CP{A(j)) = — — — ^ B+ " exp <^ -—V } + r- cxp <^ - ■? \ ■ (58) 

* norm 

aA0 v27rcrAr t 2(t^ J ^Jlixap y 2.a p J 

If we concentrate on the away-side peak, the area under the Gaussian peak Ap is proportional to the production rate 
for approximately back-to-back hadron pairs. In going from p+p to d+Au collisions, we thus have: 

~ ■ (59) 

If one knows A^p for p-t-p collisions, from the theoretically calculated in the last subsection, one can predict the 
area under the Gaussian peak A^"^" in d-|~Au collisions. 

On the other hand, the width ap oi the away-side peak can be determined from the transverse momentum imbalance 
through Eq. ()53p . The change from p+p to d+Au will mainly depend on A{q\), which can be calculated in our 
formalism from Eq. (j47p . To obtain the dihadron azimuthal correlation in d-|-Au collisions, we first fit the dihadron 
correlation in p+p collisions to obtain A^p and a^p . We get: 

B== 0.0051, Aa,= 0.0122, ctjv = 0.4238, = 0.0159, cii. 0.7135. (60) 

From Eq. using {pu_) - 2.54 GeV and {p2i_) ^ 1.28 GeV in p+p colhsion [36], we find (g^)pp - 5.0 GeV^. We 
then calculate {q\_)dAu = {qj_)pp + ^{Qa} obtain the away-side width ap in d+Au collisions. On the other hand, 
with our calculation of i?^^ we get Aj;^" from Eq. (|59p . The dihadron azimuthal correlations obtained this way are 
compared to the STAR experimental data ^] for both central and peripheral collisions in Fig. [51 The solid curves are 
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calculated with containing both dynamical shadowing and cold nuclear matter energy loss effects. The dashed 
curves are calculated with only dynamical shadowing. The constant offset is B = 0.01405 for central collisions and 
B = 0.0066 for peripheral collisions. As we can see from the plot, our calculation gives a very good description of 
the experimental data in central collisions. For peripheral collisions, the agreements get worse. The main reason for 
the deviation comes from the fact that the experimental data for peripheral d+Au collisions show a clear broadening 
effect in the away-side width ap Ji]. However, our calculated broadening A{q'j_)dAu « v4^/'^(Af^^[(6))/(A^^^i(5min.bias) 
becomes quite small. 
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FIG. 6. Azimuthal correlation associated with back-to-back dihadron production in central (top) and peripheral (bottom) 
d-l-Au collisions. Theoretical curves are calculated for (j/i) — (2/2) ~ 3.2 and (pi±) = 2.68 GeV and (p2x) = 1-31 GeV in d-|-Au 
collision Data is from STAR [61]. 
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IV. CONCLUSIONS 



In summary, by taking into account both initial- and final-state multiple parton scattering inside the nucleus, we 
calculated in perturbative QCD the increase in the transverse momentum imbalance (nuclear-induced broadening) of 
dijet and dihadron production in high energy p-l-A (d-|-A) collisions relative to the more elementary p-l-p collisions. 
The nuclear-induced broadening can be used to calculate the width of the away-side peak in dihadron correlation 
measurements. For phenomenological applications, we combined our new theoretical findings with previously derived 
coherent power correction (dynamical shadowing) and cold nuclear matter energy loss results. Perturbative QCD 
calculations that take these effects into account were recently shown to give a good description of forward rapidity 
single inclusive particle production in d-l-Au collision at RHIC. In this manuscript we provided the corresponding 
evaluation for dihadron cross sections and correlations relevant to the new STAR and PHENIX measurements. With 
cold nuclear matter parameters constrained by data on deep inelastic scattering on nuclei, we found that the calculated 
nuclear modification factor is roughly consistent with the PHENIX experimental data. Finally, by combining the 
calculated width of the away-side peak and the nuclear suppression factor, we were able to describe reasonably well 
the dihadron azimuthal correlations measured by the STAR experiment. Even though we need the baseline from p-l-p 
collisions, our formalism does describe the effects of cold nuclear matter in going from p+p to d+Au collisions pretty 
well for mid-mid, mid-forward, and forward-forward correlated hadron pairs at RHIC. 
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